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Abstract 

We present a calculation of the generalized parton distributions of the photon when there is non-zero 
momentum transfer both in the transverse and longitudinal directions. We consider only the contri- 
butions when the photon helicity is not flipped and calculate those at leading order in electromagnetic 
coupling a and zeroth order in the strong coupling a^. We keep the leading logarithmic terms as well 
as the quark mass terms in the vertex. By taking Fourier transforms of the GPDs with respect to the 
transverse and longitudinal momentum transfer, we obtain the parton distributions of the photon in 
position space. 



I. INTRODUCTION 



There has been a lot of progress recently in the study of exclusive processes like deeply virtual 
Compton scattering (DVCS) and hard meson production. The former is the process ep — )■ ejp 
where there is a momentum transfer between the initial and final proton and a real photon 
is observed in the final state At large virtuality of the inducing photon, the kinematics 
factorizes and the soft part can be written in terms of the generalized parton distributions 
(GPDs) of the proton. The GPDs are unified objects rich in content about the spatial and 
spin structure of the proton. In the forward limit when the momentum transfer is zero, GPDs 
reduce to parton distributions and their first moments give the nucleon form factors. In a 
similar manner, factorization for the process 7* ((5)7 — )■ 77 was studied in ^ for high 
and it was shown that the amplitude can be written in terms of the GPDs of the photon in 
the DVGS kinematics, that is when the momentum transfer between the real photons in the 
initial and final states is small compared to Q^. The same process in a different kinematical 
region gives information on the generalized distribution amplitudes of the photon jsj. The 
generalized distribution amplitudes are related to the photon GPDs by crossing. On the other 
hand, the parton distributions (pdfs) of the photon is of interest since a long time. They are 
anomalous as they show logarithmic scale dependence even at zeroth order in the QCD coupling 
constant ag [4j . These play an important role in photon induced high energy processes and have 
been studied extensively both in theory and experiments. Like the photon structure function, 
the photon GPDs show logarithmic scale dependence already in parton model, and are called 
anomalous GPDs. In [21] the GPDs of the photon are identified with the Fourier transform of the 
matrix elements of light-front bilocal currents between photon states, and they were calculated 
at leading order in electromagnetic coupling a and zeroth order in considering only leading 
logarithmic terms when the momentum transfer between the initial and the final photon states 
was purely in the longitudinal direction. In a previous work jsl, we calculated the photon GPDs 
when the momentum transfer was purely in the transverse direction. We expressed the photon 
GPDs in terms of the light-front wave functions of the target photon, which can be evaluated 
in perturbation theory. Fourier transform of the photon GPDs with respect to the transverse 
momentum transfer A"*- gives parton distributions of the photon in the transverse position or 
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impact parameter space. Like the proton, these have the physical interpretation of probabihty 
density [6|. In other words, the impact parameter dependent parton distribution (pdf) q{x,b^) 
of a photon is the probabihty of finding a quark with momentum fraction x at a distance b-^ 
from the 'center' of the photon in the transverse position space. Thus the photon GPDs give an 
interesting picture of the transverse shape of the photon in position space. When the skewness 
C is nonzero, there is also momentum transfer in the longitudinal direction. The photon GPDs 
obey different scale evolution equations in the regions x > ( and a; < C [2j ■ In terms of overlaps 
of the LFWFs, these can be expressed as the overlaps of two particle wave functions in the 
kinematical regions 1 > x > and —l<x<( — l. In the_ region < x < ( there is a particle 
number changing overlap, similar to the proton GPDs [7]. As shown in js] for the proton, 
when the skewness ( is non-zero, GPDs in impact parameter space do not have a probabilistic 
interpretation. However, they are still interesting as they now probe the partons when the 
initial photon is displaced from the final photon in the transverse impact parameter space. 
This relative shift does not vanish when the GPDs are integrated over x in the amplitude. In 
this work, we investigate the parton content of the photon in impact parameter space when the 
skewness is non-zero. We limit ourselves to the kinematical region 1 > x > ( and —l<x<( — l 
where only the two-particle LFWFs contribute. As shown in jsl in our previous work, photon 
GPDs in impact parameter space show distinctive features compared to the proton GPDs when 
the momentum transfer is purely in the transverse direction. In this work we show that similar 
conclusion can be drawn for a more general momentum transfer which has both longitudinal 
and transverse components. In 9f| we had introduced a boost invariant longitudinal impact 
parameter a conjugate to the skewness parameter C. It was found that the DVCS amplitude 
for a dressed electron target state shows diffraction pattern in a space. Similar pattern was 
observed in a simulated model for a meson and even in the holographic model for the meson. 
Further investigations revealed that general features of the diffraction pattern is independent 
of the model used for the meson or the proton GPDs and depend only on the finite upper limit 
of the ( integration and the invariant momentum transfer squared t. However, the appearance 
of the diffraction pattern depend on the model of the proton GPDs 10|. In this work, we 
investigate the photon GPDs in transverse as well as longitudinal position space. 

The plan of the paper is as follows : in section II we present the calculation of the photon 
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GPDs; numerical results are given in section III. Conclusions are in section IV. 



II. GPDS OF THE PHOTON FOR NON-ZERO C 

The GPDs of the photon can be expressed as the following off-forward matrix elements 
defined for the real photon (target) state j^: 



F' = J ^e^(7(P') I ^(0)7V^(y-) I 7{P))- (1) 

F'^ contributes when the photon is unpolarized and F'^ is the contribution from the polarized 
photon. We work in the light-front gauge A'^ = 0. The above GPDs were calculated when the 
momentum transferred square t = [P — P'Y was purely in the transverse direction in [5], in 
other words, when the skewness C was zero. On the other hand, in [2], these were calculated 
when the momentum transfer was purely in the longitudinal direction or A-*- = 0. We use the 
standard LP coordinates P"^ = ± P'^, = ± y^. Since the target photon is on-shell, 
P^P^ — P-*-^ = 0, the momenta of the initial and final photon are given by: 

P = ( P+ , 0^ , ) , (2) 
P' = ((1 - OP* . -A^ . (T^^) ■ P) 

The four-momentum transfer from the target is 

A = p-p' = (|cp^ A\ , (4) 

where t = A^. In addition, overall energy-momentum conservation requires A^ = P^ — P'^, 
which connects A^^, (, and t according to 

(1 - Qt = -A^'. (5) 

The details of the method of calculating P^ and F'^ using the Pock space expansion of the 
state was given in jsl in the limit of zero (. Here we do the calculation when both ( and A^ 



are non-zero. There are quark and antiquark contributions respectively for ^ < a; < 1 and 
— l<x<(^—l. In this work, we restrict ourselves to these two regions. The two-particle 
LFWFs for the photon can be calculated in perturbation theory and are given by 
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where we have used the two- component formalism ll|, ll2| and m is the mass of q{q). A is 



the helicity of the photon and Si,S2 are the helicities of the q and q respectively. The wave 
functions are expressed in terms of Jacobi momenta Xj = and q^- = — XiP^. These obey 
the relations — '^i^i^t' — 0- The GPDs can be written in terms of the overlaps of the 

LFWFs as follows : 



pi = I (fq^dxi5{x - Xi)'4j^{x, q^-il- x)A^)'^2{xi, q^) 
— j (fq^dxi6{l + X — Xi)'ip2{^j 1'^ + 5^^)^2(a^i5 q^ 



(7) 



Here x = ^3^. We have suppressed the helicity indices and the sum over them. The first 
term is the contribution from the quarks and the second is the contribution from the antiquark 
in the photon. As the light-cone momentum fraction xi has to be always greater than zero, the 
first term contributes when 1 > x > ( and the second term foT —1 < x < ( — 1. We have taken 
( to be positive. There is a particle number changing overlap in the region < x < (, Using 
the LFWFs each component can be calculated separately. We calculate in the same reference 
frame as [9|]. Note that the light cone plus momentum of the target photon is non-zero. Finally 
we get for the unpolarized photon 



F" = 5^ ^ [(1 - a;' - a; + 2x'x){h + h + Lh) + 2m^h]e{l - x)e{x - C) 



(y.s r 

^ 4^ (-a; - x" + 2(1 + + I2 + L'Q + 2m^I'^ e{l + x)e{l + x- C)(8) 



where x' = and x" = ^-jr^- Here the sum indicates sum over different quark flavors; 
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FIG. 1: (Color online) (a) Plot of unpolarized GPD F'^ vs x for fixed values of —t in GeV^ and 
at C = 0.1 and (b) polarized GPD F'i vs x for fixed values of —t in GeV'^ and at C = 0.1 and (c) 
unpolarized GPD F'^ vs C for fixed values of —t in GeV^ and at a; = 0.6 and (d) polarized GPD F'^ 
vs C for fixed values of —t in GeV^ and at x = 0.6 , A = 20GeV. 



L = -2m^ + m^x{l - x) + mV(l - x') - (A-^)^(l - x'f, L' = -2m^ - m^x{l + x) + rn^x" (1 - 
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the integrals can be written as, 
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FIG. 2: (Color online) (a) Plot of q{x,(^,b) vs x and (b) q{x,(,b) vs x for fixed values of b and at 
C = 0.1 (c) Plot of vs b and (d) q{x,C,b) vs b for fixed values of x and at = 0.1, where we 

have taken A = 20 GeV and Amax= 3GeV where Amax is the upper limit in the A integration. The 
distributions are given in GeV^ and b is in GeV 



DD' 



da- 



P(x,a,C,(AJ 



(9) 



where D = {q^f -m^x{l-x)+m'^ and D' = (g^)V(A^)^(l-x')^-2g^-A^(l-x')-mV(l- 



x')+m^, and P{x, a, (, (A-^)^) = — Q;m^x(l— x)— m^(l— x')+m^+Q;(l— x')^(A-^)^. 
As already seen in our previous paper for zero (, at zeroth order in «<, the results are scale 
dependent, this scale dependence in our approach comes from the upper limit of the transverse 
momentum integration A = Q. ju is a lower cutoff on the transverse momentum, which can be 
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FIG. 3: (Color online) (a) Plot of q{x,C,b) vs ( and (b) q{x,(,b) vs ( for fixed values of x and at 
b = 0.6 (c) Plot of q{x, C, b) vs C and (d) b) vs C for fixed values of b and at x = 0.6, where we 

have taken A = 20 GeV and Amax= 3GeV where Amax is the upper limit in the A integration. The 
distributions are given in GeV^ and b is in GeV 



taken to zero as long as the quark mass is nonzero. We included the subdominant contributions 
from the mass terms in the vertex. 

For the antiquark contributions we have similar integrals 
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FIG. 4: (Color online) (a) Plot of q{x,a,t) vs a for a fixed value of x = 0.4 and different values of —t 
in GeV^, (b) same as (a) but for x = 0.6; (c) q{x,a,t) vs a for a fixed value of x = 0.4 and different 
values of —t in GeV^, (d) same as (c) but for x = 0.6; we have taken A = 20 GeV. 
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where H = {q^)^ + m?x{l + x) + and H' = (g"*")^ + (A^ 



X 



"2 



+ 2g^ ■ AV -mV(l 



a;") + m2, and Q(x, a, C, (A-^) ) = am2x(l + a;) - (1 -a)a;"(l-a;")m2 + m2 + a(l-a)a;"2(A^) 



For polarized photon the GPD F'^ can be calculated from the terms of the form e\e] 



|2|. We consider the terms where the photon helicity is not flipped. This can be written as, 



47r2 



[x + x' - l)(/i + h + Lh) + 2rrt'h]e{l - x)e{x - Q 



Y^-l[-x- + + L'l',) + 2m'r, e{l + + X - C) (11 



In the limit A"*" = we get the expressions in [2j and the well known splitting functions 
in the forward limit. Starting from the expressions of photon GPDs, we define the parton 
distributions of the photon in impact parameter space for nonzero C, as, 



q{x, C, 6) = ^ I d'A^e-'^^-'^F'^ix, (, t). (12) 

Here 6 =| 6-*" | is the transverse impact parameter which is a measure of the transverse 
distance between the struck quark and the center of momentum of the photon. In the limit 
of zero skewness C the parton distributions in impact parameter space have a probabilistic 
interpretation 61]. However, in most DVCS experiments C, is nonzero, and it is of interest to 
investigate the GPDs in 6-*" space with non-zero As mentioned in the introduction, the 
physical interpretation of the parton distributions in impact parameter space was given for 
the proton in when the skewness C, is nonzero. In this case the transverse location of the 
target proton is different before and after the scattering. This transverse shift depends on the 
skewness C, and h and the information on the transverse shift is not washed out even if the 
GPDs are integrated over x in the DVCS amplitude. On the other hand, the boost invariant 
longitudinal impact parameter a was first introduced in and it was shown that DVCS 
amplitude shows a diffraction pattern in longitudinal impact parameter space. GPDs for the 
proton when expressed in term of a also exhibit the similar diffraction pattern |lO|. In the 
same way here we introduce the boost invariant longitudinal impact parameter conjugate to 
the longitudinal momentum transfer as a = |6^P+. The photon GPDs in longitudinal position 
space is given by: 



1 p<,max 1 f<,max 

qix, ^,t) = — dCe^'^^^'-Z'F^ix, C,t) = — dCe'<'^F'^{x, (, t) 

"1 pCtnax 1 /*Cmax 

q{x, a,t) = — dCe'<''^'-/'F^{x, C.t) = — dCe^^F^x, (, t) (13) 

2vr Jo Jq 



In the numerical calculation, we shall consider only in the region ^ < x < 1, and the upper 
limit of ( integration (max will be taken as x. 
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FIG. 5: (Color online) Plots of (a) q{x,a,t) and (b) q{x,a,t) vs a for fixed t in GeV^ and different 
values of x ; (c) F'? vs ( for — t = 3.0 GeV^ and different values of x ,(d) q{x,a,t) vs a for fixed t in 
GeV^ and different values of x for which q{x,a,t) shows two maxima with a central minima 



III. NUMERICAL RESULTS 



In all plots we have taken A =Q= 20 GeV and m=0.0033 GeV. In Figs. 1(a) and (b) we 
have plotted the GPDs for unpolarized photon F'^ and for polarized photon F'^ as functions of 
X for a fixed value of ( and different values of t. Note that in all plots x is positive and x > ( 
so that the contribution comes from the active quark in the photon. As seen in our previous 
work and also in ^ when is zero, that is when the momentum transfer is purely in the 
longitudinal direction, F'^ is symmetric with respect to a; = 1/2 when the momentum is shared 
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equally among the quark and the antiquark. This is the case when the subleading terms due 
to the quark mass at the vertex are neglected. When A-*- is non-zero, this symmetry is not 
there, the effect is more prominent at lower values of x. We see a similar effect_^when both ( 
and A"*- are non-zero. Both the GPDs and have been normalized as in [5] to compare 
with [2] in the limit of zero t. They become independent of t as x ^ 1 because in this limit all 
the momentum is carried by the quark in the photon. The slope of the polarized GPD changes 
with increase of —t and at a; ~ 0.5 it changes sign. Figs 1(c) and 1(d) show the dependence 
of F'^ and F'^ on ( for a fixed value of x and different values of t. Both and F"^ decrease 
steadily with increase of (. The rate of decrease is faster for lower values of | — t |. In Fig. 2 we 
have plotted the Fourier transform (FT) of F'^ and F^ with respect to A"*- for non-zero (. we 
took the upper limit of the A-*- integration, Amax.=3 GeV. Figs. 2(a)- 2(d) show the plots of 
the impact parameter dependent parton distributions q{x, (, h) and q{x, (, h) both as functions 
of X for fixed h as well as functions of h for fixed x; for a given value of As stated before in 
our calculation we have restricted ourselves to the kinematical region x > C,. As observed for 
C = 0, both q{x, (, b) and q{x, (, h) rise sharply as x — > 1 as the GPDs become independent of 
t. In Figs. 3(a)-(d) we have shown the behaviour of q{x,(,b) and q{x,(,b) as functions of ( 
for fixed x and and b =\ b-^ \ values. Both decrease steadily with increase of (. For a given b, 
the rate of decrease with ( increases with the increase of x and for a given x the rate decreases 
with increase of b. At a given (, the distributions are broader in b space for x ~ 1/2 when 
both q and q carry equal momenta. As stated in the introduction, the GPDs in transverse 
impact parameter space for non-zero ( probe partons inside the target photon when the initial 
photon is displaced from the final photon in the transverse impact parameter space. The x and 
( dependence of the photon GPDs play an important role in the longitudinal position space 
distribution which we show in Fig 4. Fig. 4 (a) and (b) show the plot of q{x,a,t) vs. a for 
different values of t and two different values of x. 4 (c) and (d) show the plots of g(x, a, t) vs. a. 
q{x, a, t) shows prominent diffraction pattern with a central maximum and several secondary 
maxima separated by well-defined minima at a fixed value of x = 0.4. The diffraction pattern 
is not that prominent for g(x, a, t) at the same value of x. Also for g(x, cr, t) the diffraction 
pattern becomes less prominent at higher values of x and disappears at x ~ 0.8. g(x, o", t) show 
two major peaks separated by a central minimum at x = 0.6. 
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In order to further study the diffraction pattern in a space in Figs. 5(a) and (b) we plot 
q{x,a,t) and g(x,cr, t) respectively vs a for a fixed value of t and different values of x. As 
can been seen from the plots, the positions of the first minima for both move towards smaller 
values of cr as a; increases. We compare this observation with the diffraction pattern shown by 
the proton GPDs in a space and by the DVCS amplitude for a dressed electron. There the 
positions of the first minima moved in towards smaller values of cr with the increase of —t. In 
contrast, in the case of photon GPDs the positions of the minima are independent of —t. An 
analogy with the diffraction pattern in optics was given in [ol. It was conjectured that the finite 
range of ( integration acts as a slit of finite width necessary to produce the diffraction pattern. 
(max for a given t is determined by DVCS kinematics and is proportional to (— t) for the proton 
target of finite mass. This (max is analogous to the slit width in optics experiments. There 
the width of the central maximum is inversely proportional to the slit width, and in analogy, 
in DVCS the width is inversely proportional to Cmax and therefore to —t. However the photon 
is massless and as seen from Eq. ([5]), for a photon target (max is 1 in the limiting sense when 
(A^) = for a given t. In our analysis, (max = x. So the slit width in the optics analogy is 
proportional to x and as expected the minima move to smaller values of cr as x increases. As 
seen in Figs. 4(d), 5(b) and 5(d), q{x,a,t) has two maxima separated by a central minimum 
for those values of x for which F'^ changes sign in the ( range covered as seen in Fig. 5 (c). At 
X = 0.6, F"^ changes sign almost in the middle of the ( range covered : it is positive for half 
of the ( range and negative for the other half as a result q{x, a, t) has the lowest value for the 
minimum at x = 0.6. Also as seen from the figs. 5 (c) and 5 (d), the central minimum for q is 
not observed when F"^ is either positive or negative during most of the ( region covered. Thus 
the pattern in a space depends closely on the x, ( dependence of the photon GPDs. The height 
of the maxima of both q{x, a, t) and q{x, a, t) decrease with increasing —t for a given x. 



IV. CONCLUSION 



In this work, we have calculated the GPDs of the photon when the momentum transfer 
has both transverse and longitudinal components, in other words, when both ( and are 
non-zero. We calculated both the polarized and unpolarized GPDs in terms of the light-front 
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wave functions of the photon which, in turn, can be calculated in perturbation theory. The 
helicity of the photon may flip in the process. In the present work we considered only the 
non-helicity flip terms. We calculated them at leading order in electromagnetic coupling a and 
zeroth order in strong coupling a^. We also kept the subleading mass terms in the vertex. The 
GPDs depend logarithmically on the scale. However, here we studied the x, C and t dependence 
at a fixed scale. By taking a Fourier transform with respect to A-*- we expressed the GPDs in 
the transverse impact parameter space. As already observed for zero C, the photon GPDs show 
distinctive features compared to the proton GPDs in impact parameter space. For non-zero C 
similar featured are observed. We introduced a boost invariant longitudinal impact parameter 
a conjugate to Taking a FT with respect to C we expressed the GPDs in longitudinal position 
space. They show an interesting pattern similar to the diffraction pattern in optics. This was 
also observed for the GPDs of the proton as well as in the DVCS amplitude for a dressed 
electron. We presented a comparative analysis of the behaviour of the photon GPDs and that 
of a proton in a space and showed that the finite range of the C integration as well as the x, C 
dependence of the photon GPDs determine the pattern in the longitudinal position space. As 
is well known, the substructure of the photon is of relevance in high energy processes. Here we 
showed that the photon GPDs can provide a unique picture of the structure of the photon in 
3D position space. Further investigation would involve a study of the photon GPDs when the 
helicity is flipped. It would be interesting to understand if the GPDs of the photon obey the 
general properties like polynomiality and positivity. A violation of positivity was seen in the 
transverse impact parameter space in [sl and that needs to be investigated further. One has to 
remember also that a complete picture of the GPDs in a space would need the contribution from 
the region x < The analysis presented here provides a major step towards understanding 
these objects. 
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